Nonreciprocal transmission of electromagnetic waves by three-layer magneto-optical mediums
Yu Guan-Xia1, 2, †, Fu Jing-Jing2, Du Wen-Wen2, Lv Yi-Hang2, Luo Min1
College of Science, Nanjing Forestry University, Nanjing 210037, China
College of of Information Science and Technology, Nanjing Forestry University, Nanjing 210037, China

 

† Corresponding author. E-mail: sys@njfu.edu.cn

Abstract
Abstract

We investigate the non-reciprocal transmission properties of a three-layer structure filled with magneto–optical medium and normal medium. Based on the transfer matrix method, we deduce the total transmission coefficient for a one-dimensional (1D) structure with anisotropic mediums. When two-side layers with magneto–optical medium loaded in opposite external magnetic field, the time-reversal symmetry of transmission properties will be broken. Our numerical results show that the non-reciprocal transmission properties are influenced by external magnetic fields, incident angle, and thickness of the normal medium layer. Since the non-reciprocal properties can be easily realized and adjusted by the simple structure, such a design has potential applications in integrated circulators and isolators.

1. Introduction

The phenomena of magneto–optic effect (MOE) discovered by Faraday in the middle of the 19th century,[1] which triggered a flood of research on fundamental physics and applications due to its special properties,[27] such as the Faraday effect, Cotton Mouton effect, and magneto–optical Kerr effect. Because the electromagnetic (EM) parameters of magneto–optical mediums are anisotropic and asymmetric, and also can be modulated by the external magnetic field, one of the other novel properties in the magneto–optical mediums (MOMs) is nonreciprocity.[813] In other words, the wave propagation in the MOMs breaks the time-reversal symmetry. Because of its important applications in integrated circulators and isolators, the nonreciprocal properties of MOMs have attracted considerable attention.

However, the weakness of MOE leads to the obstacles of magneto–optical effects in practical applications. Therefore, it is essential to sweep the obstacles to seek for the strategy to enhance MOE. Surface electromagnetic waves, which exist at the interfaces separating different media, are proven to be important approaches in enhancing the MOE. One of the most well-known types of these waves is optical Tamm states (OTSs) which have existed in a photonic crystal.[1417] At the interface between two photonic crystals with overlapping bandgaps, there are surface electromagnetic OTSs that exhibit nonreciprocal behavior in magneto-photonic structures. Furthermore, a similar confined electromagnetic modes can also exist at the boundary between a metal layer and a Bragg mirror, which are called Tamm plasmon polaritons (TPPs).[18] When the medium in the Bragg mirror is a gyrotropic material, TPPs have nonreciprocal properties. The other types of such waves are surface plasmon-polaritons (SPPs),[19,20] which exist at interfaces separating metals and dielectrics. When the metal is MOMs, the time-reversal symmetry of the SPs is broken and the nonreciprocal wave propagation can be observed.[21,22] Although these two types of surface waves have different natures, their properties are quite similar.

In the magneto-photonic crystals, nonreciprocality of propagating waves can be switched by the one-way electromagnetic OTSs (or TTPs), which is based on an infinite or half infinite periodic structure. In reality, the properties of nonreciprocality could be diminished by the finite periodic structure. At the same time, the nanoscale magneto–optical layers are tightly linked in the magneto-photonic crystals, it is extremely difficult to magnetize magneto–optical layers by the external magnetic field. Furthermore, the anisotropic EM parameters in the magneto–optical mediums are not constant, and are regulated by the magnitude of magnetic field and frequencies of incident waves. The nonreciprocality is intensively related to EM parameters in the magneto–optical mediums, which have been seldom involved in previous researches. Therefore, it is necessary to search effective methods to realize and enhance nonreciprocality.

In this paper, a simple three-layer structure composed of anisotropic MOMs magnetized in the opposite directions and normal mediums is constructed, and the nonreciprocal propagation property is enhanced by nonreciprocal SPPs. The rest of this paper is organized as follows. First, the three-layer structure with MOMs is constructed, and reflection coefficient is conducted by the transfer matrix method.[23,24] Second, the numerical results and discussions associated with our purposes are presented. Finally, a conclusion is given in the last section.

2. Physical model and computational method

A one-dimensional (1D) structure with three layers is illustrated in Fig. 1, where A and layers are anisotropic MOMs, and C layer is the normal material silicon dioxide (SiO2). In the structure, dA, , dC are the thicknesses of layers A, , and C, respectively. B represents the external magnetic fields on the MOMs, where dot and cross represent the magnetic field direction along the positive and negative y directions, respectively. In the Voigt geometry, the relative permittivity and permeability of layer A and can be expressed as

Layers A and C are magnetized plasma materials, which the relative permittivity can be expressed as
where ω is the incident wave angular frequency, and is the oscillation frequency of electric plasma, which is related to the number density of electrons n, the electronic charge e, electronic mass m, and the permittivity in the vacuum. is the gyrofrequency, which is related to the magnitude and direction of the additional static magnetic field B. If the direction of additional magnetic field is along the y direction, then we take the minus sign in Eq. (2). Conversely, we take the plus sign in the formula.

Fig. 1. The schematic diagram of 1D structure air/ /air, composed of plasma layers A, , and normal material layer C.

For simplicity, a TM mode incident plane wave, , incident on the stratified medium, the total field in one layer can be written as:

and represent the wave propagates along the positive and negative z directions, respectively, and their fields are related to those in neighboring regions by boundary conditions. According to the transfer matrix methods,[23,24] the fields of boundary in one region i is connected to those in the adjacent region j by the transfer matrix:
where Tij is the transfer matrix of the adjacent layers i and j, can be written as
for anisotropic medium, where , , , and
is the wave vector of z axis in the layer l, and kx is the transverse wave vector along x direction. If the medium is isotropic, Nl and Ml can be simplified by and . When the TM mode wave propagates through this multilayer structure, the incident, reflected, and transmitted electric fields are connected via a total transfer matrix
where Pi is the phase-shift factor of the one layer from front boundary to rear boundary, and defined as follows:
where di is thickness of the i-th layer. The total transmission and reflection coefficients t and r of the EM field through such a structure is given by
where M11 and M21 are elements of the transfer matrix M.

According to Eq. (2), the Voigt permittivity in the gyroelectric medium, the function of the frequencies, is less than zero when the frequencies of incident waves are less than certain frequencies. Under this circumstance, the behaviors of anisotropic gyroelectric medium and SPPs in the boundaries between gyroelectric medium and normal medium can be stimulated. Furthermore, the asymmetry of the gyro-parameters in two-sides gyroelectric medium is enlarged by the opposite magnetic fields, which enhance the nonreciprocal transmission of the electromagnetic waves.

3. Numerical results

To demonstrate the nonreciprocity in the designed structure, the transmission properties are numerically calculated when a plane wave is incident from left boundary of the structure. The medium parameters in the numerical simulations are chosen as: For medium in layer C, we choose material SiO2, the relative permittivity ϵr = 2.07, the relative permeability μr = 1; for gyroelectric medium in layers A ( ), n = 8 × 1017 m−3, , , and relative permeability μr = 1. Because the Voigt permittivity of gyroelectric medium is less then 0 when the frequencies of incident waves are less than 8 GHz, in the simulations we chose to use a frequency range from 2 GHz to 8 GHz.

3.1. The influence of the incident angle on nonreciprocity

First, we assume that a TM wave is incident on the PCs structure at different angles θ = 60°, 45°, and 30° from the boundaries, and we choose the external static magnetic field B = 0.02 T, the thickness of , dC = 30 mm. Figure 2 shows the relationships between transmission coefficient and frequencies, where the solid line and dash line correspond to the forward ( ) and backward ( ) incidences, respectively. For incident angle θ = 60° in Fig. 2(a), there are two different transmission peaks, f = 5.10673 GHz for positive incidence and f = 5.22093 GHz for negative incidence. The results show that there are transmission nonreciprocity due to the asymmetric structure A and in the opposite external magnetic fields. Figure 2(b) and 2(c) also show two different transmission peaks for incident angle 45° and 30°. As the incident angles decrease, there is also transmission nonreciprocity, and the two different transmission peaks shift right and broaden.

Fig. 2. The relationships between the transmission coefficient and the frequencies for different incident angles, where the solid line represents positive incidence, and dot line represents negative incidence. (a) incident angle θ = 60°, (b) incident angle θ = 45°, and (c) incident angle θ = 30°.

To clearly show the transmission nonreciprocity, the field patterns of Hz for the forward ( ) and backward ( ) incidences are shown in Fig. 3, where the incident angle is 60°, and other conditions are the same as given in Fig. 2. Figure 3(a) and 3(b) show the field patterns of Hz for frequency of f = 5.10673 GHz, when counter-propagating TM plane waves are incident from left and right boundaries. For the positive incidence, multi-layers look like a transparent medium, and the fields can completely pass through them. At the same time, there are strong fields in the interface between gyroelectric medium and normal medium and in the middle normal medium. In contrast, for the negative incidence, the transmitted waves are almost suppressed, and the fields cannot transmit through the mutli-layer. Then, the stationary waves in the incident region are formed by interactions between incident waves and reflective waves. For the frequency of f = 5.22093 GHz, the similar results are shown in Figs. 3(c) and 3(d), the negative propagating wave can transmit through the multi-layers and the positive incident waves are completely reflected. Thus, the nonreciprocal propagating properties are realized by the simple three layers.

Fig. 3. The distributions of magnetic fields Hy along xz plane, where the external magnetic fields B = 0.02 T and incident angle is 60°. (a) Positive incidence with f = 5.10673 GHz, (b) negative incidence with f = 5.10673 GHz, (c) negative incidence with f = 5.22093 GHz, and (d) positive incidence with f = 5.22093 GHz.

To further investigate the nonreciprocal transmission properties in the multi-layers, the distributions of magnetic fields Hy along the z direction with x = 0 has performed for incident angle 60° in Fig. 4, where the solid line and dot line represent and incidences, respectively, and the dash line represents the different layers. For positive incident frequency f = 5.10673 GHz, when the waves transmit through the muti-layers, as figure 4(a) shows: there are strong SPPs in the interface between the gyroelectric medium and normal medium, and there appears resonant peak in the middle normal medium. For negative incidence at the same frequency, there are no SPPs in the interface and no resonant peak in the normal medium, which appears to be a one-way total transmission property due to nonreciprocity of the gyroelectric medium. For frequency f = 5.22093 GHz, the resonant peak appears in the multi-layers for negative incidence and the positive incidence is almost completely suppressed due to the nonreciprocity of gyroelectric medium.

Fig. 4. The distributions of magnetic fields Hy along the z direction with x = 0, where the external magnetic fields B = 0.02 T and incident angle is 60°. (a) f = 5.10673 GHz, the solid and dot lines represent positive incidence and negative incidence, respectively, (c) f = 5.22093 GHz, the solid and dot lines represent negative incidence and positive incidence, respectively.

In the three-layer structure , the gyroelectric medium A and on both sides of normal medium C constitute two asymmetric units AC and . At our working frequencies, the Voigt permittivities of the gyroelectric medium are negative, and the waves cannot be transmitted in the propagating direction due to imaginary propagation constant. But in a certain frequency, the unidirectional SPPs (also known as surface magneto-plasmons (SMPs)[22]) in the surface between gyroelectric medium and normal medium can be excited by the external fields in the bandgap, then the incident wave can completely propagate through the three-layer structure. On the contrary, when the wave at the same frequency is incident from the opposite direction, the SPPs are completely compressed or even removed in the interface, and incident wave is completed reflected.

3.1.1. The influence of magnetic field B on nonreciprocity

Because the nonreciprocity of the structure is related to the external magnetic fields, we have investigated the transmitting properties of multi-layers for different magnetic fields, where the incident angle is 60°, and the other conditions are the same as given in Fig. 2. When the external fields B increase from 0.2 T to 0.6 T, the gyroelectric components of gyroelectric medium largen, and nonreciprocity of structure is enhanced. The relationships between the transmission coefficient and the frequencies for different external magnetic fields are depicted in Fig. 5, where the solid line and dot line represent and incidences, respectively. This clearly shows the transmitting peak shifts to low frequencies for case, and the transmitting peak shifts to high frequencies for case, when the external magnetic field increases. That is to say, the distance of frequencies of counter-propagating waves are symmetrically broaden, and the nonreciprocal transmission properties increases with the increase of the external magnetic fields.

Fig. 5. The relationships between the transmission coefficients and the frequencies for different external magnetic fields with the incident angle θ = 60°. The solid and dot lines represent positive and negative incidences, respectively. (a) B = 0.02 T, (b) B = 0.04 T, and (c) B = 0.06 T.
3.2. The influence of the thickness of layers on nonreciprocity

To take into account the thickness of layer C between two gyroelectric mediums, we calculate the transmission coefficients and the results are shown in Fig. 6 for the case of different thicknesses of layer C. Based on this discussion, the behaviors of two gyroelectric mediums are similar to those of metal, and a resonant cavity is constructed. The different thicknesses of layer C decide the different resonant frequencies. Figure 6 shows the the relationship between transmission coefficients and frequencies for TM mode plane wave vectors and . When the thickness of layer C decreases from 30 mm to 20 mm, for positive incidence, the transmission peaks shift from f = 5.10673 GHz to f = 6.60446 GHz; for the negative incidence, the transmission peaks shift from f = 5.22093 GHz to f = 6.81596 GHz, and the distance of transmission peaks for counter-propagating waves broadens, which means that the nonreciprocity of the structure is improved with the decrease of thickness for layer C. Therefore, if the layer C is consisted of multilayer with the same thickness and EM parameters, then a specific nonreciprocal transmission frequencies can be acquired by increasing or decreasing layers, which have potential applications in integrated circulators and isolators. In addition, because the layer C is a normal medium, in the hamburger structure , the properties transmission nonreciprocity can be realized by layer C filled with other EM parameters.

Fig. 6. The relationships between the transmission coefficients and the frequencies with different thicknesses of layer C, where the solid and dot lines represent positive and negative incidences, respectively. The incident angle θ is 60°, and the external magnetic fields B is 0.02 T. (a) dC = 30 mm, (b) dC = 25 mm, and (c) dC = 20 mm.
4. Conclusion

In this paper, the transmission nonreciprocity has been studied by a simple structure that consists of three layers filled with MOM and normal medium. When two layers with magneto–optical medium loaded in opposite external magnetic field, the time-reversal symmetry has been broken. Numerical results show that the nonreciprocal transmission peaks are influenced by the external magnetic fields, incident angles, and thickness of the normal layer. With the increasing of external magnetic fields and incident angle, the nonreciprocal transmission peaks shift to high frequencies, and the properties of the transmission nonreciprocity are enhanced. At the same time, with the decrease of thickness of the normal medium, the nonreciprocal transmission peaks also shift to high frequencies and the properties of the transmission nonreciprocity are also enhanced. Because the nonreciprocal properties can be easily realized and adjusted by the simple structure, this design has potential applications in integrated circulators and isolators.

Reference
[1] Faraday M 2005 Experimental Researches in Electricity 3 614
[2] Yablonovitch E 1987 Phys. Rev. Lett. 58 2059
[3] John S 1987 Phys. Rev. Lett. 58 2486
[4] John S Wang J 1990 Phys. Rev. Lett. 64 2418
[5] Shiveshwari L Mahto P 2016 Solid State Commun. 138 160
[6] Wu C J Chen D X King T C Lin W C 2015 Appl. Opt. 54 8738
[7] Wang Y Zhang D G Xu S X Xu B G Huang T 2017 Micro. Opt. Technol. Lett. 59 1347
[8] Yu Z F Veronis G Wang Z Fan S H 2008 Phys. Rew. Lett. 100 0203902
[9] Fu J X Lian J Liu R J Gan L Lia Z Y 2011 Appl. Phys. Lett. 98 211104
[10] Shoji Y Mizumoto T 2014 Sci. Technol. Adv. Mater. 15 506
[11] Leviyev A Stein B Galfsky T Krishnamoorthy H Kuskovsky I L Menon V Khanikaev A B 2017 Appl. Photon. 2 3966
[12] Christofi A Stefanou N 2018 Phys. Rev. 97 125129
[13] Wu Z Chen J G Ji M X Huang Q Z Xia J S Wu Y Wang Y 2015 Appl. Phy. Lett. 107 221102
[14] Khanikaev A B Baryshev A V Inoue M Kivshar Y S 2009 Appl. Phy. Lett. 95 011101
[15] Dong H Y Wang J Cui T J 2013 Phys. Rev. 87 045406
[16] Fang Y T Chen L K Zheng J Zhou L Y Zhou J 2014 IEEE Photon. J. 6 4801611
[17] He C Sun X C Zhang Z Yuan C S Lu M H Chen Y F Sun C 2013 Opt. Express. 21 28933
[18] Kaliteevski M Iorsh I Brand S Abram R A Chamberlain J M Kavokin A V Shelykh I A 2009 Phys. Rev. 76 165415
[19] Economou E N 1968 Phys. Rev. 182 539
[20] Li J F Liang H B Wang L 2017 Chin. Phys. 26 114103
[21] Yu Z Veronis G Wang Z Fan S 2008 Phys. Rev. Lett. 100 023902
[22] Shen L F You Y Wang Z Y Deng X H 2008 Opt. Express 23 950
[23] Fang Y T Ni Y X He H Q Hu J X 2014 Opt. Commun. 320 99
[24] Fang Y T 2014 Plasmonics 9 1133